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We present a statistical readout method for quantum interferences based on time series analysis of 
consecutive single electron transfers through a double quantum dot Aharonov-Bohm interferometer. 
Waiting time distributions qualitatively indicate the presence of interferences and provide informa- 
tion on orbital-detuning and coherent interdot-electron transfer. Interdot transfer induced oscilla- 
tions are Aharonov-Bohm phase sensitive, while those due to level detuning are phase-independent. 
The signature of the quantum interference in the waiting time distribution is more apparent for 
weakly coupled electron transfer detectors. 
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Double-quantum dot (DQD) junctions provide an ex- 
perimental setup to study phase coherent transport 
[Hi 0, 9 ^ii'i to realize Aharonov-Bohm (AB) interfer- 
ometers 0, H 0] • This is of general interest as they are 
potential candidates for qubits. A crucial aspect for their 
realization is the noninvasive determination of the pres- 
ence of quantum interferences in order to minimize deco- 
herence. 

So far theoretical studies on transport properties of 
double quantum dot AB interferometers have been fo- 
cused on average current^ i, 1, [13, [O, [H, [11 [3 and 
shot noise [1^, [la, [l3, [3 properties. Recently time- 
resolved detection of single electron transfers in single 
quantum dots [l^, [20, l2l| and DQD in series has become 
experimentally feasible ^^- Waiting time distributions 
of consecutive electron transfers can be obtained from 
time-series analysis and provide detailed information on 
quantum dots [23, [2^ and single molecules [25]. They 
were found to be sensitive to interference due to mul- 
tiple electron paths in DQD junctions [23] and contain 
more detailed information than current and noise mea- 
surements |24| . 

In this letter we propose a weakly invasive statistical 
method based on waiting time distributions of single ex- 
ternal electron transfers that can determine the presence 
of quantum interferences, small detunings of the DQD 
orbitals, inter-dot transfer coupling, and Coulomb inter- 
action. These quantities are connected with qualitatively 
distinguishable oscillations in the waiting time distribu- 
tion. These oscillations are sensitive to the AB phase, 
(j) = $/<&o, and are suppressed bX (p — rnr for an integer 
n when interdot transfer is present. Here, $ is the mag- 
netic flux perpendicular to the junction and $0 = h/e 
the magnetic flux quanta. In contrast, the oscillations 
purely due to energy detunings of the DQD orbitals are 
(^-independent. We show that their detection requires 



weakly coupled electron detectors thereby avoiding to 
inflict fast decoherence on the DQD which qualifies the 
proposed method as a readout scheme for coherently op- 
erating qubits. 

To that end, we exploit a master equation in the many- 
body Fock space of the DQD, assuming weak system- 
reservoir coupling. For simplicity we consider spinless 
electrons and each QD dot can hold only one electron 
at most. We decompose the total Hamiltonian of the 
DQD-AB interferometer junction into Ht = Hs + Hr + 
Hsr{4>)- The DQD part (system) reads 
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A(clc2 + 4ci)- (1) 



Here, e^ with s = 1 or 2 is the orbital energy of the 
specified QD; U denotes the strength of the Coulomb 
repulsion between electrons; A is the inter-dot electron 
transfer parameter in the DQD-AB interferometer. The 
Hamiltonian of the electrodes is given by two indepen- 
dent free electron reservoirs, Hr = X],y=i,r Z]g ^qvc\vCqy 
The index v denotes the left [l) or right (r) electrode, 
q their intrinsic degrees of freedom. The electron cre- 
ation (annihilation) operators c\ and cjj^ (c^ and Cg^) 

satisfy the anticommutator relations: {cfe,cj,, } = Skk' 
and {c[.,cj,,} — {ck,Cki} = 0, for all k,k' = s,qv. The 
system-reservoirs coupling responsible for electron trans- 
fer between the electrodes and the DQD reads 



^sk(0) = E E [t^':H<p)cW+^.c. 



(2) 
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The AB phase-dependent transfer parameters satisfy 
T^l\<j,) = r,ie^*-/2 and T^l\^) = T^2e-''^"'^ for the 
two parallel dots pierced by a magnetic flux considered 
here, with (pi = —(p^ = </< to account for the sign change 
in the phase between coupling to left and right electrode. 
We describe the DQD by the reduced density operator 
p{t) of the system. Note that the inter-dot electron trans- 
fer A in Eq. ^ leads to off-diagonal elements in the sys- 
tem Hamiltonian Hs in the orbital basis. We transform it 



into eigenbasis, Hs = 0~^HsO, where O consists of the 
eigenvectors of Hg. The same transformation is apphed 
to the creation and annihilation operators ^^ = 0~^CsO 
and *t ^ 0-140, so that Hsr = O'^HsrO. The stan- 
dard perturbation theory leads to the quantum master 
equation |23l |: 



pit) = -iCp{t) + J2 (-n. + s+ + s;)p(t). 



(3) 
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The system Liouvillian £ ■ — [Hs , ■ ] describes the coher- 
ent dynamics. The dissipative superoperator in Eq. ([3]) is 
separated into the diagonal contribution 11^ that leaves 
the number of electrons in the system unchanged, and the 
off-diagonal E+ and S~ for the increase and decrease the 
number of electrons in the DQD, respectively; see Ref.[23l 
for the derivation. This separation is necessary in order 
to keep track of the trajectories of single electron trans- 
fers. Consider, for example, the scenario of detecting an 
electron entering the DQD through the left electrode at 
time to and leaving through the right electrode at time 
t. The waiting-time distribution of consecutive electron 
transfer events is then given by the joint-probability [23|: 



Pi^r{t,ta)^ti-s{J:;St. 



,T,+ps{to)}, 



(4) 



with St.to = exp[{—i£ — Hi — Ilr) (t — to)] being the 
propagator of the system in the absence of transfer events 
at the electrodes within the waiting time interval. Quan- 
tity (HI) can be obtained from the time-series of single 
directionally resolved electron transfers between the elec- 
trodes and the system. One has to record a sufficiently 
large number of the I — > r events and generate a his- 
togram of the number of occurrences as function of the 
time interval t — to- The histogram has to be normalized 
by the total number of considered events. 

The aforementioned physically distinct dissipative 
components are formally given as 23']: T,'^ = 

and n^ = J2s (^l^i-t^ + ijiil^ + H.C.). The involved 
superoperators are defined as the left- or right -actions 

(vjf . = v]/ . or vl> ■ = ■ vp) of the associated Hilbert-space 
^-operators. Besides the annihilation (creation) opera- 
tors_3[c (vl'l) we also have to consider their auxiliaries 

lillii: 



KVit,cb) 



dt'Cl'^lit-t';q^)e-'^^'-''^^s'. (5) 



Here, d+J, (t; cj>) = E,,> T^s^*i^)Tl'lm^UtWAO))n 

and C^-J,{t;^) = T,,'T^'^\<t>)TP; {<j^){c,.{t)cl,M) H 
are the AB phase-dependent interacting reservoir cor- 
relation functions. Applying the given phase rela- 
tions in Tqs {(f)) and assuming further Tqi — Tq2 — 
Tq for the AB phase-free parts lead to the relations: 
Cl'^ht)^^ Citkt)^ = Ci-^\t), &^l{t) = d^\t)e^*^, 
and Cl,2i{t) = Gv ' {t)e~^'^" . The auxiliary operators in 



their non-Markovian form [Eq. ([5])] can be numerically 
evaluated numerically without further approximations as 
shown in Ref . [23 and [26l 

To derive analytical results we apply the Born-Markov 
approximation, together with the wide-band limit for the 
reservoir spectral density. The latter leads to Cl (t) — 
r/;°d6#)(6)eT-*. Here, /+(6) = 1 - /-(e) = 
[e(/3(e-/^,.) _^ i]-i = j(g _ ^^) is the Fermi distribution 
function, with /3 being the inverse temperature and [in 
the Fermi energy of the electrode v. The Born-Markov 
approximation amounts to replacing the range of time in- 
tegration in Eq. ([S]) with (—00,00). As results, the auxil- 
iary annihilation operators defined in Eq. (0) can be eval- 
uated as 



*lt^ = r/±(/:)(*i + *2e 
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(6a) 
(6b) 



which depend on the AB-phase {(j) ~ (f>i — ~(j)r) but no 

long on the time. The auxiliary creation operators '^Is 
are of similar expressions, but with the replacements of 
£ -^ -£, -^ -(/) and ^, -^ *| in Eq. ©. Note that 
since Hs is diagonal in the many-body Fock space, the 
action of the superoperator /^(£), which is determined 
by the Fermi function and the diagonal system Liouvil- 
lian, can be carried out easily. All the 16 auxiliary op- 
erators, ^Is and ^li with v — l,r and s = 1,2, can 
now be evaluated [cf. Eq. ^] in terms of 4x4 matrices in 
the Fock-space representation. Consequently, the action 
of each dissipative tensor in the second term of Eq. ([3]), 
which has been given in terms of the left- and right- 
multiplications of some \l/s (\1/J) and ^i-s {^i/s ) is now 
determined. It is worth to mention here that the approx- 
imation scheme explored in Eq. ^ leads to an Eq. ^ in 
Lindblad form. 

We use the following parameter scheme to describe our 
calculation results. A bias of 2V is applied symmetrically 
/i;/r = Mcq ± V . The Orbital energies of the DQD are set 
to be ei = Eg+a and £2 = Eg ^ ct; i-e., the orbital energy 
split (or detuning) is 2a. We set the vacuum DQD state 
eo = as the energy zero, and tg — 1 the internal energy 
unit. In all calculations, /icq = 1.0 and T — 0.1. 

Figure [1] demonstrates the dependence of the wait- 
ing time distributions Pi-,r{t) (left-panels) on the AB- 
phases (j), together with their Fourier transforms F(uj) 
(right-panels) exemplified at three representing values of 
4> — 0, 7r/2, and tt. The Coulomb repulsion parameter 
U — 1.0 and the bias 2V = 0.2 are common, while the 
interdot transfer and orbital energy split parameters are 
(A, a) — (0.1, 0) in the upper, (0, 0.05) in the middle, and 
(0.1, 0.05) in the bottom panels, respectively. Clearly, the 
influences of A and a on the waiting time distribution 
are qualitatively distinct, especially in the two limiting 
regimes. While Pi^r{t) shows only little dependence on 
(j) in the dot orbital-split case (the middle panel: A = 
but a 7^ 0), it is strikingly sensitive to the AB-phase in 
the interdot-transfer case (the upper panel: A 7^ but 
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FIG. 1: Pi^r{t) (left-panels) as functions of AB-phase (j) 
and time t; The corresponding Fourier transformation F{uj) 
(right-panels) a.t (j) ~ (dash), 7r/2 (solid), and n (dot), re- 
spectively. The upper, middle and bottom panels are for three 
representing sets of interdot-transfer rate A and orbital de- 
tuning a. Other parameters are U — 1.0, T — 0.1,2V = 0.2 
and jj.cg = 1.0 (in unit of eg); see text for details. 



a — 0). In the latter case, the characteristic oscillation 
is maximized at = 7r/2, but disappears at (f> = and 
(j) = TT. These observations can be largely understood as 
follows. 

The dot orbital-split (A = but a ^ 0) case resembles 
the transport through double slits. The resulting inter- 
ference [23| persists and is insensitive to the AB phase, 
due to the fact that (j)i = ~4'r in each orbital channel 
largely cancels out the AB-phase effect. This accounts 
for the basic feature observed in the middle panels of 
Fig.ffl 

In interdot-transfer case (upper panels: A 7^ but 
a = 0), the aforementioned double-slit feature is de- 
stroyed. The interdot-transfer allows electrons to switch 
between the two pathways provided by the DQD. Thus, 
different phases can be accumulated as the electron trans- 
fer through the coupled DQD and the aforementioned 
phase symmetry is broken along some of the possible 
transfer trajectories. As a result, the total accumulated 
phase depends on the value of 0. It leads to a pure decay 
of Pi^r{t) at the AB phase = or tt, where e*"^ = e"*"^. 
However, at other values of (j>, it leads to an effective 
phase difference between the eigen-levels which are sub- 
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FIG. 2: Pi^r{i) as function of inter-dot transfer rate A and 
time t and corresponding Fourier transformation F{uj). The 
AB-phase is <j) = 7r/2 and the orbital detuning a = 0. Other 
parameters are same as Fig. [1] 



ject to an induced energy gap of 2A, responsible for the 
AB-phase activated oscillations observed in the upper 
panels of Fig.[TJ 

In the intermediate regime shown in the bottom panels 
of Fig.[l] oscillations can be observed for all 0; however, 
the Fourier transform reveals a frequency shift when AB- 
phase is tuned by the magnetic field. At (/) = 7r/2, the 
observed frequency corresponds to the DQD eigenenergy 
gap (2VA2 + a^ = 0.224), while at </> = or tt, it is blue 
or red shifted, respectively. The amplification of the os- 
cillation at = 7r/2 is characteristic for inter-dot transfer 
and allows to distinguish it from orbital detuning. The 
latter causes only small oscillations at = or tt. 

Figure[2]examines further the influence of A on Pi^r{t) 
(left) and its spectrum F{oj) (right), with a — and 
(j) = 7r/2, where oscillations due to AB phase-activated 
interferences between the eigen-levels are at maximum. 
Note that the interference would remain dark at </> = in 
this case; as can be seen the upper panels of Fig.[TJ The 
amplitude of Pi^r{t) oscillation decreases with A, which 
corresponds to a decreased average current through the 
DQD. 

To analyze other coherent operation conditions, let us 
focus on the orbital-detuning (A = and a 7^ 0) case 
where Hs is diagonal. We also neglect the Liouville- 
space off-diagonal elements in 11; -I- 11^, which have a 
relatively small influence in the weak coupling regime. 
As results, the propagator St^to in determining Pi^r{t, io) 
[Eq. (|3])] becomes diagonal, and the analytical solution is 
achievable. Moreover, the waiting time distribution is 
separable into Pi^r{t) = Posc{t) + Pdecay{t)- A detailed 
discussion of the decaying terms Pdecay (t) that depends 
weakly on the AB-phase would exceed the scope of this 
letter; but it has been provided for the case of incoherent 
transport through single benzene molecules [25l| which 
can be applied to QD-systems as well. As the coherent 
operation conditions are concerned, we focus only on the 
oscillation term, which is independent of the AB-phase 
for the orbital-split case (cf. the middle panels of Fig. [1]). 

Poscit) = Por^b^iV, a)e-2«(c/,v')rt^Qg(2at), (7) 
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FIG. 3: The damping parameter a{U, V) and the pre- 
exponential amphtude parameter b {V,a), as their function 
dependence, for the oscillation term Posc{t) [Eq. ([7])]- Inter- 
dot electron transfer is absent A = 0. 



where po is the initial vacuum state occupation number, 
a{U, V) = a{U, -V) ^ f{U + V) + f{U -V) + 2, and 



KV,a) 



1 + e^"^ + 2e(^+°)^ 

g/3(y+2a+l)/2 



f{V + a)f{V-a). (8) 



The left-panel of Fig. [3] depicts the damping parameter a 
as function of U and V. It assumes the maximum value of 
4 for small Coulomb coupling U < V and is independent 
of a. Also the decay rate proportional to the system- 
electrode coupling strength F. Thus a weak coupling is 
required for the observability of interferences. This qual- 
ifies statistical analysis of waiting time distributions as 
an indirect method to study internal processes indirectly 
avoiding fast decoherence in the system. 

The right -panel of Fig. [3] depicts the pre-exponent pa- 
rameter b^ as function of V and a. It reveals further 
the parameter regimes where oscillations are observable. 
One condition is that V < a. Oscillations are suppressed 
at negative bias larger than the DQD energy gap. The 
amplitude is strongly increased when V > a. However in 
this regime the decay rate 2aF may reach its maximum 
and prevent the observability of coherence. Apparently, 



the presence of strong Coulomb coupling, as well as oper- 
ating at small bias regime, are favored for the observation 
of interference effects by means of waiting time distribu- 
tions. 

In conclusion, a Markovian quantum master equation 
in the Fock space was formulated and employed to calcu- 
late the waiting time distribution of consecutive electron 
transfers in AB interferometers. Based on this we de- 
scribe a novel statistical method to determine quantum 
interferences, inter-dot electron transfers, orbital detun- 
ing and the AB-phase. Orbital detuning and inter-dot 
transfer induce oscillations in the waiting time distribu- 
tion in the presence of interference. The two cases can 
be distinguished qualitatively since the latter one is sensi- 
tive to the AB-phase. The observability of oscillations re- 
quires the presence of strong Coulomb interaction, small 
bias and a weak electrode-system coupling. 

The indirectness of the statistical detection avoids fast 
decoherence but a large number of transfer events is nec- 
essary in order to extract information. This might be 
advantageous for a qubit in operation with a continu- 
ous readout. The method does not provide information 
on a single operation but can determine whether a large 
set of operations is carried out coherently. Also other 
sources of decoherence like coupling to phonon bath have 
to be minimized. The signature of interferences in wait- 
ing time distribution can survive in the presence of a 
phonon-bath [2J|. The proposed scheme can be realized 
utilizing presently available technology. For that pur- 
pose two DQD in series junctions which act as detectors 
by their coupling to their respective quantum point con- 
tacts [22| should be installed on both sides of a parallel 
DQD junction. This setup consisting of six QDs avoids 
decoherence inflicted by the charge state measurement of 
the quantum point contact. 
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